Geometrodynamic Quantization 
and Time Evolution in Quantum Gravity 



Arkady Kheyfets* 
Department of Mathematics, 
North Carolina State University, 
Raleigh, NC 27695-8205 

Warner A. Miller 
Theoretical Astrophysics Group (T-6, MS B288), 
Theoretical Division, 
Los Alamos National Laboratory, 
Los Alamos, NM 87545 
(February 7, 2008) 

We advance here a new gravity quantization procedure that explicitly utilizes York's analysis 
of the geometrodynamic degrees of freedom. This geometrodynamic procedure of quantization is 
based on a separation of the true dynamic variables from the embedding parameters and a distinctly 
different treatment of these two kinds of variables. While the dynamic variables are quantized 
following the standard quantum mechanical and quantum field theoretic procedures, the embedding 
parameters are determined by the "classical" constraint equations in which the expectation values of 
the dynamic variables are substituted in place of their classical values. This self-consistent procedure 
of quantization leads to a linear Schrodinger equation augmented by nonlinear "classical" constraints 
and supplies a natural description of time evolution in quantum geometrodynamics. In particular, 
the procedure sheds new light on the "problems of time" in quantum gravity. 
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I. INTRODUCTION. 



The task of describing time evolution in quantum gravity appears to encounter numerous and seemingly insur- 
mountable obstacles. The difficulties are so pervasive that some researchers have voiced a doubt whether the concept 
of time and time evolution can be adequately introduced in quantum gravity. The issue of time evolution in quan- 
tum gravity has been reviewed by Kuchaf [g| who formulated it as a set of "problems of time." He concluded that 
the problems of time have not been successfully resolved in either the Dirac or the square-root ADM approaches to 
gravity quantization. In this paper we review and analyze both these procedures of quantization in light of a new 
approach to gravity quantization. |^ 

The classical theory of gravity is a fully constrained theory. It is our thesis that this feature should not be mirrored 
in quantum gravity, and our observation that the source of the "time" difficulties in both the Dirac and ADM 
quantization approaches arises because the York split [Q of the dynamical variables from the time parameterization 
variables occurs too late. Here we propose to introduce such a split before the dynamical theory is developed. This, 
we suggest, will supply a more cogent quantization procedure as well as bypass the "problems of time." 

We provide in (Sec. II) a brief review of the standard techniques of canonical gravity quantization. In Sec. Ill 
we advance a new procedure of gravity quantization (geometrodynamic quantization). This procedure incorporates 
York's analysis of the geometrodynamic degrees of freedom and appears to avoid the difficulties of the standard 
procedures of quantization. In Sec. IV we provide two simple examples to illustrate our geometrodynamic procedure of 
quantization; namely, we quantize the axisymmetric Kasner and Taub cosmologies. Sec. V provides a brief discussion of 
the specific features of this quantization approach, while Sees. VI and VII provide a brief description of the "problems 
of time" in standard canonical quantum gravity and the solution of these problems within our new framework. 
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II. CANONICAL QUANTUM GRAVITY. 



Both the Dirac and square-root Hamihonian approaches to canonical quantization make use of the 3+1 spht of 
the spacetime geometry by Arnowitt, Deser and Misner (hereafter referred to as ADM). According to the ADM 
prescription of classical general relativity one considers a slicing of the spacetime by a family of spacelike hypersurfaces 
labeled by a parameter t. This parameter can be thought of as a time coordinate so that any slice is identified by 
the relation t = const. The remaining three spacetime coordinates determine a coordinatization of each slice. The 
spacetime metric {^'^^g^.u) is parameterized by the shift (A^*), lapse (N) and the 3-metric of the slice (gik)- 
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The standard action of general relativity for the gravity field 



(2.1) 



after its augmentation by appropriate boundary terms can be represented as 



(2.2) 



(2.3) 



In this equation tt'-' , 7i, and TC' are given as follows 
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where K^^ is the extrinsic curvature tensor, K = Kl is the trace of the extrinsic curvature tensor, and 11 is the matrix 
of TT*^ . All the quantities in (2.4) - ( |2.6| ), including the covariant derivative in (2.6) are related to the 3-geometry of 
the slice. 

Without any further refinement or redirection, the ADM approach treats all six components of the 3-metric gik 
as the gravitational "field coordinates" (with their conjugate momenta tt*''), while the lapse N and the shift 
are treated as Lagrange multipliers (no conjugate momenta for them are present in the expression for action). The 
expression 



nADM{^'\g^3,N,N,) = Nn{TT'\g,,) + N,n\TT'\g,,) 



(2.7) 



plays the part of Hamiltonian when variations of tt*^ and gij are considered. These variations produce the twelve 
Hamilton equations. 



dt 
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dHADM 



(2.9) 



We do not display this last equation in all detail as it is rather lengthy and is not necessary for this paper. In addition, 
there are three supermomentum constraints 
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(2.10) 



and the super-Hamiltonian constraint 



H(7r'^5y) = 0. 



(2.11) 



These four con st rain t equ ati ons a re obtained by variations of Ni and N in the ADM action, respecti vely . The 
equations (2^), (^.9D, (2.9), ( 2.10 ) are not independent. The interplay between the Hamilton equations (|2.8|), ( |2.9| ) 
and the constraints ( ^.IC ), (2.11) is rather involved. In particular, if the constraints are satisfied on an initial slice 
then the Hamilton equations guarantee that they are also satisfied on all spacelike slices. On the other hand, if 
the constraints are satisfied on each slice of all possible spacelike foliations of a given spacetime then necessarily 
the Hamilton equations are satisfied. This last feature of general relativity has led people to refer to gravitational 
dynamics as a "fully-constrained theory." While it is clear that this property can be valuable in the classical theory, 
it is also equally clear that it is implicitly based on the assumption of the uniqueness of the spacetime 4-geometry 
This last assumption of uniqueness becomes problematic in quantum geometrodynamics, and consequently, the 
fully-constrained property of geometrodynamics also becomes problematic. This situation is, in our opinion, the root 
of the notorious "problems of time" in quantum gravity. 

A further investigatio n of t he ADM picture of gravitational dynamics leads one to the interpretation of the super- 
momentum constraints (2.1C) as expressing the 3-dimensional diffeomorphism invariance (the freedom of choice of 
coordinates on the slices) and the super-Hamiltonian constraint (2.11), together with the natural identification M 



5S 



(2.12) 



as the Hamilton-Jacobi equation of the theory 



(2J3) 



Although this Hamilton- Jacobi equation involves all six components of the 3-metric, the 3-diffeomorphic invariance 
as expressed by the supermomentum constraints. 



5S 



0, 



(2.14) 



allows one to identify the Hamilton-Jaco bi eq uation as an equation that describes the evolution of the 3-geometry 
rather than of the 3-metric. Equation ( 2.13| ) has numerous disturbing features. The most disturbing of them is 
that, according to it, the super-Hamiltonian Ti cannot be interpreted as the generator of the time translation 
In addition, unli ke th e standard situation in mechanics, the super-Hamiltonian Ti. participating in the Hamilton- 
Jacobi equation ( 2.13 ) does n ot co incide with the TLadm Hamiltonian (2/7) which generates the dynamic evolution 



via Hamilton equations ( 2 



This Hadm Hamiltonian is related to the Lagrangian C by the standard relation 



Ti-ADM — TT 



dt 



C 



(2.15) 



while, ordinarily, TL is not. 

Dirac's procedure of canonical quantization utilizes the standard prescription 



5S 

Sgij 
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^ij = -7 

I dgij 



(2.16) 



and leads to four functional differential equations based solely on the four constraint equations of the classical theory 



W {n'^,g,,)^ = 0, 



(2.17) 



> 9ij 



(2.18) 



The three equations (2.17) are often interpreted as a requireme nt for the state functional to be a functional of the 
3-geometry — ^![^^'Q\ rather than of the 3-metric. Equation ( 2. 18] ) is the Wheeler-DeWitt equation [0,0 which 
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is considered to be a proper wave equation for quantum gravity. This equation reminds one more of a Klein-Gordon 
equation than a Schrodinger equation. As we have already mentioned, the state functional 5* in this equation, 
although originally introduced as a functional of the 3-metric 5" = \E'[5'i/c], after imposing on it the requirement of 
3-dimensional diffeomorphic invariance, is considered to be a functional of the underlying 3-geometry 5" = ^'[^^•'fj]. 
The commutation relations are imposed on all the components of the 3-metric, although proper care is taken to 
ensure that they are diffeomorphically invariant. The Wheeler-DeWitt equation, just as its classical counterpart, the 
Hamilton-Jacobi equation, does not admit an interpretation of the Hamiltonian as a generator of time translations. 

An attempt to remedy this situation via introducing, instead of the super-Hamiltonian TL a Hamiltonian that, in a 
sense, does generate the time translation, has been undertaken in the ADM square-root procedure of quantization . 
In this procedure the set of six 3-metric parameters is split in two subsets /32}, and {ai, a2, as, ri} of parameters. 
The first of them is treated as the set of true gravitational degrees of freedom, and the second is consid er ed as the 
set of embedding variables. The classical super-momentum and super-Hamiltonian constraints ( [2.10 ), (2.11) are 



then solved on the classical level with respect to the momenta conjugate to the embedding variables, which, after 
natural identifications (for the sake of simplifications in notations we omit indices on (3 and a parameters; they can 
be recovered easily whenever it is necessary) [[zljiol] 



leads to the following four equations: 

g = -,.,(|,Aa,n) (2^20) 

(2^21, 

One of the functions /iq, contains a square root operation in it (we assume from now on that it is h^), which is 
why it is called the square-root Hamiltonian. The standard prescription for quantization, 

T 'i^i^, T 'i^i^, T ^Tf^^-^, (2.22) 

Sa 6a sn Sn' (5/3 i Sf3' ^ ^ 

leads to the four quantum equations, 

ih^ = ha(n'^,l3,a,n)'^, (2.23) 
da 

in^^^hn{^^,P,a,n)M>, (2.24) 

where = ^ [/3,Q;,r2]. The choice of parameters a and Q is frequently made in such a way t hat the three equations 
(2.23) are interpreted as a requirement of the diffeomorphism invariance, while equation ( ^.24 ) is considered to be the 



proper Schrodinger equation. In this case parameters a can be interoreted as coordinatization parameters while the 
function n can be thought of as the many-fingered time parameter |ll| , p^ with hn treated as the Hamiltonian. 

We wish to point out that the square-root Hamiltonian hn in the classical theory, just as it was the case with the 
super-Hamiltonian 7i, is not related to the Lagrangian C as it would be in a standard dynamic theory (cf. equation 



Both the ADM square-root and Dirac quantization procedures admit, in principle, a time evolution prescription to 
be introduced in them, although this should be done on different stages. In Dirac's quantization procedure the state 
functional ^I* is a functional of either the slice 3-metric or the slice 3-geometry, which can be written in terms of /3, a, 
n functions as 

* = *P',g] /3] or * = ^-P^] = /3], (2.25) 

where a, n and f3 are treated on equal footing. In particular, a computation of any quantity similar to an expectation 
function, if it were possible, ordinarily would involve functional integration over either all six functional parameters a, 
ri, /3 or the three functional parameters f3 and Along these lines we wish to remind the reader here that, in Dirac's 
quantization approach, the commutation relations are imposed on all the components of the 3-metric, and in quantum 
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mechanics the evaluation of expectations involves integration over all the variables participating in the commutation 
relations. Such a calculation, however, would ordinarily exclude a possibility for one to introduce a functional time 
and a time evolution in the theory. A way to avoid such an unpleasant situation is to split the set of variables in 
two subsets, one of which (/?) is to be considered as a set of true functional arguments of the state functional (and 
to be integrated over during computing the expectations), while the others (a, Q) are to be considered as a set of 
functional parameters (one per slice) and to be interpreted as a representation of functional time. This split should 
be introduced only after the solution for the state functional is obtained and prior to the evaluation of expectation 
values. Such a procedure is thoroughly artificial and amounts to treating the embedding variables as quantum ones 
in one part of the theory yet as classical ones in another. 

The square-root Hamiltonian quantization is based on the same classical picture of the evolution in geometrody- 
namics as a change of 3-metric or 3-geometry from one spacelike slice to another and leads one to a Hamilton-Jacobi 
equation. However, it introduces the split of variables just before quantization. From this moment on it assumes 
that there is only one function fl per slice. The state functional 5* when restricted to one slice becomes essentially a 
functional of /? only, which can be properly written as 

^^^[a,n;(3] or ^^^[r?;/?]. (2.26) 

The computation of expectation functions (or any other momenta) now assumes functional integration over the two 
parameters /? only, with the result depending on as a parameter (of course, it will also depend on three other 
functional parameters a, which we call coordinatization parameters). This fimctional parameter (Q) plays the role 
of the many-fingered time parameter. To summarize, in the square-root Hamiltonian approach the switch from 
one picture of evolution to another occurs at the final stage of formulating classical geometrodynamics, and comes, 
just as in Dirac's approach, as an afterthought aimed to meet the difficulties generated by the original picture of 
evolution (as evolution of the 3-geometry) in the quantized theory. Such a switch does not cause any problems in the 
classical theory, but comes at a price in quantum geometrodynamics by generating numerous difficulties, including 

the notorious "problems of time" (cf. Sec. VI). 

The equations of quantum gravity, be it equations ( ^.17 ), ( 2.18| ) of Dirac's quantization or equations (2.23), (2.24) 



of the square-root Hamiltonian quantization, do not provide for a natural split. The split is, to an extent, arbitrary. It 
reflects, for any particular gravitational system, our understanding of the system's dynamics. In a sense, it is similar 
to the situation described by N. Bohr regarding the split of any considered phenomenon into its quantum and classical 
parts. Such a split depends not only on the object that we are considering but, also, on the questions that we are 
asking. 

It appears to us that the equations of quantum gravity themselves do not imply the necessity of a picture of temporal 
evolution for quantum gravitational systems. The assumption of a temporal evolution should be added, at some level, 
so as to not contradict to the rest of the theory. The first step in this direction is to introduce a split of metric variables 
into the true dynamic variables P on one hand, and the coordinatization parameters a and the many-fingered time 
parameter Q on the other hand. The variables (3 are the true quantum variables (calculation of expectation values 
involves functional integration of 5* only over these functions), while a and Q are merely functional parameters that 
allow one, after some manipulations, to introduce the concept of a temporal evolution. The two different approaches 
described in this section introduce the split of variables at different stages. This leads to difficulties of both the 
technical and conceptual nature in quantum geometrodynamics. 

In the next section we advance an alternative procedure of quantization - a procedure that, from the very beginning, 
is based on a picture of geometrodynamic evolution induced by York's analysis of the geometrodynamic degrees of 
freedom. This alternative approach does not involve a change in the paradigm of time evolution in geometrodynamics. 



III. GEOMETRODYNAMIC QUANTIZATION. 



We propose here an alternative approach to the quantization of gravity. Our approach is based on the post- 
ADM achievements made in classical geometrodynamics. jl^ In particular, we are referring to York's solution of the 
initial-value problem and his analysis of the gravitational degrees of freedom. This development was initially 

motivated by Wheeler's semi-intuitive remark that the 3-geometry of a spacelike hypersurface has encoded within it 
the two gravitational degrees of freedom as well as its temporal location within spacetime. It is this notion that the 
3-geometry is a carrier of information on time that has been referred to as "Wheeler's many-Bngered time." |]lT| , |l2| ,[7| 
It was J. York who first made this thesis precise. He forwarded what has now become almost the canonical split of the 
3-geometry into its underlying conformal equivalence class (its shape representing the two dynamic degrees of freedom 
of the gravitational field coordinate per space point) and the conformal scale factor (its scaie representing Wheeler's 
many-fingered time). Only the conformal 3-geometry is truly dynamic in the sense that it can be specified freely as 
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the initial data. The scale factor is non-dynamic and essentially specifies Wheeler's many-fingered time. The results 
of York have demonstrated that the true dynamic part of the gravitational field is not the 3-geometry but only its 
conformal part, and that the proper configuration space or "arena for geometrodynamics" should be the underlying 
conformal superspace (the space of all conformal 3-geometries) rather than Wheeler's superspace (the space of all 
3-geometries) . The conformal scale factor and three other functional parameters of the 3-metric (responsible for 
coordinate conditions) thus become external parameters. In what follows we associate the true dynamic variables (3 
with conformal 3-geometry, the many fingered time variable U, with the 3-geometry scale factor, and the remaining 
three variables a of the 3-metric with a coordinatization of a spacelike 3-surface. The many-fingered time variable 
and coordinatization variables are assumed to be fixed initially by some conditions, and then subsequently controlled 
by the shift and lapse. 

Our proposed procedure of gravity quantization is based, from the very beginning, on York's analysis of gravitational 
degrees of freedom. We suggest that one should interpret geometrodynamics as an evolution of the conformal 3- 
geometry (3 in an external field determined by the scale factor f2 and coordinatization variables a. Such an approach 
calls for a reformulation of geometrodynamics on the classical level. We start from the standard Lagrangian C (written 
in terms of the 3~metric, shift and lapse) and the associated action (with appropriate boundary terms, as needed, to 
remove the second time derivatives terms) and we introduce the momenta conjugate to the true dynamic variables 

We then use these tt/j's to form the geometrodynamic Hamiltonian TLdyn of our approach, 

Ti-DYN = T^pl3 ~ C. (3.2) 

The new Hamiltonian JioYN is distinctly different from TLadm and its arguments do not coincide with those of 
T-Ladm, namely 

Ti-DYN ='HDYN{^,Oi;P,Tri3)- (3.3) 

The variables preceding the semicolon are treated as describing an external field, while the ones following the semicolon 
are the coordinates and momenta of the gravitational true degrees of freedom, i.e. of the conformal geometrodynamics. 
The variation of /3 and tt^ leads to the equations of geometrodynamics, i.e. to two pairs of Hamilton equations, 

P-^^, (3.4) 



dUp 

dp 



YN 



(3.5) 



and, subsequently, to the Hamilton-Jacobi equation 



5S_ 

5t 



-Ti-DYN U^, a; /? 



5(3 



Here 5 is a functional of (3 and, in addition, depends on the same parameters as Ti-oYN, 

S^S[n, a; 13] . 



(3.6) 



(3.7) 



Neither the Hamilton equations (3.4), (3.5) nor the Hamilton-Jacobi equation ( |3.6| ) are capable of providing any 
predictions as their solutions depend on the functional parameters fl and a which are not yet known. One can 
complete the system of equations by adding to the Hamilton equations, or to the Hamilton-Jacobi equation, the 
standard constraint equations of general relativity. They should be satisfied when the solution for (3, of equations 
of conformal geometrodynamics (with appropriate initial data) is substituted in them (we use symbols [/3]s, [T^fj\s for 
such a solution) 



W {n,a,[l3]sA-^p]s)^Q 



(3.8) 
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The resulting equations are equivalent to the standard equations of classical geometrodynamics. It should be empha- 
sized that these constraint equations cannot be derived from variational principles in our theory. 



For the purpose of quantization, we start from our Hamilton-Jacobi equation (3.6), describing effectively what 
we refer to as conformal geometrodynamics. That is the evolution of the dynamic variables corresponding to the 
conformal part of the 3-geometry, where this evolution is parameterized by the "external" field represented by (1) the 
scale parameter and (2) the coordinatization variables. Using the Hamilton-Jacobi equation ( |3.6| ) we may transition 
to the corresponding Schrodinger equation 

in— = nDYN{.^,a-p,TTp)^ (3.9) 

where Tip = jjp- The Schrodinger equation ( j3.9| ) treats the scale parameter and coordinatization functions as external 
classical fields and quantizes only the true dynamic variables, p. The state functional ^> in this equation is a functional 
of (3 and also depends on the functional parameters and a. 

* = * a; /3] (3.10) 

This Schrodinger equation (with specific initial data) can be solved (cf., for instance the example of the Bianchi lA 
cosmological model in Sec.^. The resulting solution of this Schrodinger equation is not capable of providing any 
definite predictions as it depends on four functional parameters fi, a which remain at this stage undetermined. All 
expectations, such as, for instance the expectation values of /? 



< [3 {-^sm^s) = J ^*J^sVp (3.11) 

or of Tff3 

< TTp >,= {^s\nfj\^s) - J "^Inf^^sVP (3.12) 

also depend on these functional parameters. To specify these functions we resort to the constraint equations. As in 
case of classical geometrodynamics the constraints should be imposed on the solution of the initial-value problem of 
conformal geometrodynamics, and in this way, determine the unique values of fl and a. It is possible that there are 
several ways to couple the constraints to the quantization of the true dynamic variables, /3. We propose here that the 
four constraints be imposed only on the expectation values of the conformal dynamics 



W{^,a,<P>s,<Tr(3 >.) = 

n{n,a,< (3 >s,< >s) = 



(3.13) 



i.e. only on measurable quantities. In so doing, we explicitly avoid the interpretational conundrums associated with 
the problems of time, and we form a "classical" gravitational clock driven by the quantized gcomctrodynamic system 
- i.e. quantum-driven many-Gngered time. 

This procedure of quantization utilizes explicitly the correct treatment of the geometrodynamic degrees of freedom 
and introduces a meaningful time parameterization utilizing the shift and lapse. The super-momentum and super- 
Hamiltonian constraints are not satisfied as strictly as in the Dirac and ADM square-root approaches. Our point of 
view is that the question of slicing independence of evolution is not a well posed question in quantum gravity and 
should be recovered only in the classical limit. It is clear that this weakening of the constraints, so that they hold 
only on the expectation values over the solution of the initial-value problem of the Schrodinger equation, makes the 
theory less restrictive and enlarges the set of possible solutions. 

The goal of this paper is limited to a clarification of our basic thesis regarding the quantization of gravity — 
quantum geometrodynamics. We believe that the best way to achieve this goal is to consider a couple of simple 
illustrative examples. 



IV. QUANTUM GEOMETRODYNAMICS: THE KASNER AND TAUB COSMOLOGIES 

A. Quantum Geometrodynamics of the Bianchi lA Cosmology. 

The Bianchi 1 A cosmological model is commonly referred to as the axisymmetric Kasner model | |l5| . Its metric is 
determined by two parameters, the scale factor H. and the anisotropy parameter f3 (we choose TV* = and = 1 
values of shift and lapse for this example) 
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As this cosmology is homogeneous the two functions VL and (3 are the functions of the time parameter t only. The 
scalar 4-curvature can be expressed in terms of these two functions to yield the Hilbert action and, after subtracting 
the boundary term, the cosmological action, 

*/ 

Ic^Ih + ^^^^-'%'o =^J{p'- ^''""dt, (4.2) 

to 

where V = J J J dxdydz is the spatial volume element. 

We treat the scale factor Q,{t) as the many-fingered time parameter and the anisotropy (3{t) as the dynamic degree 
of freedom. The momentum conjugate to j3 is 

P,= ^-^ = ?^e--/3. (4.3) 

(as it is usually the case with homogeneous cosmologies, we are working here with the momentum pp rather than 
with the density tt^). The Hamiltonian of the system in our approach can be expressed in terms of the momentum 
conjugate to /3 and the Lagrangian. 

Hdyn = PpP - L 

= |^e3"f^| + ^f2^e-3^. (4.4) 

In the classical theory this Hamiltonian can be used to produce either one pair of Hamilton equations or the equivalent 
Hamilton-Jacobi equation. In any case, the dynamics picture derived in this way is incomplete. To complete it we 
impose the super-Hamiltonian constraint. 

^2 _ f^^\ „-6f2fS2 

Using the Hamilton-Jacobi equation, 



-g^=-HnYN[^g^Mt),m], (4.6) 



together with the expression (4.4) for the Hamiltonian H and the standard quantization prescription we obtain the 



Schrodinger equation for the axisymmetric Kasner model 

The constant h in this equation should be understood as the square of Planck's length scale, rather than the standard 
Planck constant. We wish to stress here that the scale factor f2 in the Schrodinger equation is so far an unknown func- 
tion of time. This means that the equation does not describe completely the quantum dynamics of the axisymmetric 
Kasner model. To complete the dynamics picture we follow our prescription and impose, in addition to equation ( [I.7| ), 
the super-Hamiltonian constraint. 

<P/5>?=(^)'c-^"fi^. (4.8) 

Here < p/3 >s is the expectation value of the momentum = j-^ 

00 

<pp>,^(^4pf,\^SJs)^ I ^:iP,t)pp^4P,t)dp (4.9) 
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where is the solution of the Schrodinger equation with specified initial data. The system of equations (4.7), (4.S) 
provide us with a complete quantum dynamic picture of the axisymmetric Kasner model and, when augmented by 
appropriate initial and boundary conditions, can be solved analytically. 

Before discussing the initial value conditions we will find the general solution of the Schrodinger equation considering 
the scale factor as a function of time generating an external potential. For this we separate variables. 



After substituting ( 4.1C ) in the Schrodinger equation (4/7) we obtain, 



ifupT 



,30 



3V 



T4)" + —ife-^^^Tct). 



where the prime means differentiation with respect to [3. Rewriting it as 



= —ihc 



-30 



T W 



-6no2 



= -A, 



3V (I) T 8tt 

where A is the constant of separation, we obtain the equations for (t>{(3) and T[t). 



2Trh 



rA0 = 



(4.10) 



(4.11) 



(4.12) 



(4.13) 



Equation ( 4.13| ) admits only positive eigenvalues for A. Introducing the notation = fc^ we can write the solutions 

.^fc(/3), rfe(t) for fee (-00,00). 

n{t) = exp |- J f (^k^ + ^e-e^fl^") e^^dtX (4.15) 



Using the superposition of these solutions we come up with the general solution of the Schrodinger equation (4.7). 

*(/3,t)= / A,e^^^^e^v[-\j^ (^|^fc2 ^ l^e'^^li^^ e^^dtj dfc (4.16) 

— 00 

To specify a particular problem one has to furnish appropriate initial data. 

00 

*(/3,0|to =*(/3,io) = j Ake^^^dk (4.17) 

— 00 

It can be done either by specifying a function ^o) ^^id then recovering Aj^ from the equation 

00 

*(/3,to)= / Aue^'^l'dk, (4.18) 



using Fourier transforms, or by assigning Ak as a function of fc, depending on the type of the problem to be formulated. 
In this section we consider the simplest example comparable with the quantum mechanics of a particle, namely a 
wave packet. To describe a Gaussian wave packet centered initially at the value fcg of ^ (we will describe the meaning 
of fco later) we assign 

Ak ^ Ce''''^^-^"^\ (4.19) 

where the constant a effectively determines the initial width of the wave packet in momenta and C is the normalization 
constant. This leads to the following expression for the initial values of the wave function: 
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oo 



The value of the normahzation constant C is determined by the condition 

a . 



(4.20) 



(4.21) 



which leads to the value of 



Using expression (4.19) for Ak and introducing notations for f and g, 



/ = fit) 



27r 
W 



(4.22) 



3V 



n^e~^^dt, 



we can write down the solution \E's(/3,i) for the wave packet. 

oo 

*.(ft<, = ae-*./e-*-«-e*»e*'-... 

— CSO 

After a simple transformation this expression can be rewritten in the form, 

*,(/?, i)=Cexp<j-[(/3-fco/)fco-.g]S' / e-'^'=%-^-^'='e^('3-2'=''-^)'=&. 



(4.23) 



(4.24) 



(4.25) 



The integral on the right hand side of (4.25) can be evaluated. The final expression for the solution describing a 
Gaussian wave packet may be written in the following form which will prove to be convenient for future calculations: 



2 , P 



a {f3-2koff 
4 ( a2 + |; 



exp S t(/3 - kof)ko > ■ exp < i 



(/3-2fco/)2 



where, 



cos(26') = a/Ja^ + p /h\ sin(26') = {f/K)/do? + p/Ti 



exp <, ~j^g ~ It 



F2 /t;2 



(4.26) 



(4.27) 



Although expression ( 4.25| ) looks quite involved the last three exponential factors are phase factors and do not 
complicate the determination of the expectation values of the observables. 

It is clear that this solution of the Schrodinger equation describing the wave packet cannot provide any definite 
predictions as it contains the two functions of time f{t) and g(t) which are themselves related to the as yet undeter- 
mined scale factor fl. To find Q{t) we need to (1) compute the expectation < Pf} >s of the momentum pfj = j-^, (2) 
substitute this expectation value into the constraint (4^) and (3) solve the resulting equation with respect to We 
start from computing < p/s >«. 
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exp ' 



dp = ko. 



(4.28) 



In other words the expectation value of the momentum < pp >s does not change with time. It is determined by the 
/c-center of the packet at f = to- Substitution of this resuh in (6.5) yields 



P - 
ftp — 



3V 

An 



(4.29) 



This equation and the classical equations are identical. Therefore, we need not desc ribe it in detail. We only wish 
to point out once more that after the solution of this equation is substituted in (4.26) the geometrodynamic problem 
(4.7), ( |4.8D for the wave packet (4.19) is solved completely. To summarize, the many-fingered time of quantum 
geometrodynamics in case of a Gaussian wave packet of axisymmetric Kasner spacetimes coincides with its classical 
counterpart if the expectation value of the momentum of the packet is identified with the (conserved) value of the 
momentum of the classical solution. 

The expectation value for the anisotropy parameter /3, where /3 is the only quantum dynamic variable in this model, 
is given by: 



P exp 



if3-2kof)' 



d/3 = 2fco/(t). 



(4.30) 



Thus "the center" of the wave packet evolves as the classical Kasner universe determined by the momentum value 
equal to fco would evolve. The spread of the wave packet with time is the variance in /3. 



<{p~<p>r >, 



iP - 2koff exp . 



{P-'^koff 



dp = 



(4.31) 



It is obvious from ( [4.3l| ) that the spread of the packet increases with time. The result is similar to that of the quantum 
mechanics of a free particle; after all the Bianchi I cosmology is the free-particle analogue of quantum cosmology. 



B. Quantum Geometrodynamics of the Taub Cosmology. 



We present within this section a second application of our approach to quantum gravity. In addition to illustrating 
our theory on a relatively simple model, we show that very concept of time emerges by imposing the principle of 
general covariance as weakly as possible. In particular, the Hamiltonian constraint is imposed as an expectation- 
value equation over the true dynamic degree of freedom of the Taub cosmology - a representation of the underlying 
anisotropy of the 3-space. In this way the concept of time appears to be inextricably intertwined and woven to the 
initial conditions as well as to the quantum dynamics over the space of all conformal 3-geometries. This quantum 
geometrodynamic approach will ordinarily lead to quantitatively different predictions than either the Dirac or ADM 
quantizations, and in addition, our approach appears to avoid the interpretational conundrums associated with the 
"problems of time." Q So without further ado, in this subsection we apply our quantum geometrodynamic approach 
to the Taub cosmology jist and numerically solve the coupled Schrodinger and expectation- value constraint equations. 

The Taub cosmology is an axisymmetric homogeneous cosmology parameterized by a scale factor il(t), and an 
anisotropy parameter P{t). The line element may be expressed as 

ds^ = -dt^ + al^"" {e'P)^^ aV^ (4.32) 

where (/3) = diag{P, P,—2P), and the one forms as, = cosipdO + sinipdip, = sinipdd — costpdcf), and cr"^ — 
dij} -f cos 9d(j). The scalar 4-curvature is expressed in terms of $7 and P and yields the action, 

37ro3 r ■ 1 (3) 

Ic^^J W -^')-^ R}e'"dt. (4.33) 
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Here ^^^R 



(4 — e represents the scalar 3-curvature. 



We treat here the scale factor D,{t) as the many-fingered time parameter |ll|] and the anisotropy P{t) as the dynamic 
degree of freedom. W The momentum conjugate to f3 is obtained from from the Lagrangian, C 



Pf) 



dC_3 
9^ ~ 2 

The dynamical Hamiltonian for this cosmology can be expressed in terms of this momentum and the Lagrangian 



(4.34) 



1 



l(3) 



R 



(4.35) 



In the classical theory Hdyn can be used to construct either the Hamilton- Jacobi equation or the two Hamilton 
equations; however, to complete the dynamics we must in addition impose the super-Hamiltonian constraint, 



l(3) 



R 



(4.36) 



Using this dynamical Hamiltonian (Eq. ( [4.35 )), the corresponding Hamilton- Jacobi equation and the standard 
quantization prescription we obtain the Schrodinger equation for the Taub cosmology. 



ih 



di 



2^ 'djP 



(4.37) 



The scale factor Q.vaV and m in this equation should be treated as an unknown function of time, t. To complete the 
quantization we impose, in addition, Eq. (4.36) as an expectation-value equation o ver [3 (where < • >s= J ^* • ^ sd(3 
and s is the solution of of ( 4.37 ) augmented by appropriate initial data, cf. Eq. ( 4.39 )). 



\ m J 



]:<^'^R>s 
6 



(4.38) 



This system of equations Eq. ( 4.37 ) and Eq. ( 4.38 ) provide us with a complete quantum-dynamic picture of the 
axisymmetric Taub model, and when augmented by the appropriate boundary conditions can be solved numerically. 
We display in Figs. 1-2 a solution of the these two equations for an initially Gaussian wave packet. 



(4.39) 



with = y277r/A/3, po = -100, ?i = 1, /3o = 0, A/3 = 0.1 and 17(t = 0) = 1. The solution was obtained using a 
2nd-order split operator unitary integrator which preserves the norm exactly (up to round off errors) . u& 
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FIG. 1. The quantum geometrodynamics of the Taub cosmology. The solid lines in each of the four graphs represent the 
quantum solution, while the small circular dots in the upper two graphs represent the classical trajectory. The scale factor's 

(n) dependence on time {t) is shown in the upper loft, while the expectation of the anisotropy (< /3 >) is displayed in the 
upper right. The lower left shows the variance in ^ s throughout its evolution, and the lower right graph is a snapshot of the 
potential iV) at the end of the simulation. 
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FIG. 2. The Taub wave function. A snapshot of the wave function at three distinct times throughout the evolution is 
represented on each of the three graphs. The first row plots ^'^^'s, while the second and third are the real and imaginary 
parts of ^'s, respectively. The right-most wave function (plotted in a slightly narrower line width) shows the initial [t — 0) 
gaussian wave function. The left-most and higher-peaked curve represents the wave function near the bounce ai t — 1, while 
the remaining curve displays the right-going wave function well after the bounce at t = 4. This refiected wave function is 
spreading in time and therefore appears, and will continue to appear, lower in amplitude than the initial wave function. 

We have introduced in this section an apphcation of our approach to the quantization of gravity (this approach 
need not be restricted to quantum cosmology). When the Taub cosmology was quantized our theory generated a 
quantum picture based on a post- ADM treatment |Q of the gravitational field dynamics and appears to be free from 
the conceptual difhculties related with time usually associated with both the Dirac and ADM quantization procedures. 
Only the dynamical part, the underlying anisotropy (/3) has been quantized. The Hamiltonian participating in the 
quantization {Hdyn) is not a square-root Hamiltonian. This absence of a square-root Hamiltonian is a generic feature 
of our quantum geometrodynamics procedure. 

We have demonstrated here using the Taub cosmology that the concept of time may very well be inextricably 
intertwined and woven to the initial conditions as well as to the quantum dynamics over the space of all conformal 
3-geometries. 
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V. THE CHARACTERISTIC FEATURES OF THE QUANTUM GEOMETRODYNAMIC PROCEDURE. 

The quantum dynamical picture described in the two previous sections has two features that are not encountered 
typically in most of the more common quantum dynamical schemes. They should be kept in mind in order to avoid 
errors and misinterpretations in applying our procedure. 

The first such feature is related to the structure of the Hamiltonian. Formally, the Hamiltonian appears as a 
Hamiltonian of a system placed in external time-dependent field, at least when the Schrodinger equation is analyzed. 
However, the external field is determined via the constraints by the quantum state of the system. Ordinarily, the 
Hamiltonian depends on the initial data. This feature is not unique for our approach as it is encountered elsewhere 
in standard quantum mechanics of Hartree-Fock systems. 

The second feature is a split of the system parameters in two groups: (1) the truly dynamic variables (to be 
quantized); and (2) the descriptors of a "natural observer" associated to the system, driven as they arc by the quantum 
geometrodynamics. Such a split is introduced prior to quantization, and is crucial to our theory. This feature can 
be clearly observed in both the general description of the geometrodynamic quantization procedure (Sec. Ill) and in 
the examples of Sec. IV. Essentially, this split is a part of our solution of the time problems. Most of the difficulties 
related to the problems of time are caused by an attempt to enforce general covariance at the quantum level. Such 
an attempt for gravitational systems is bound to fail. It leads to the questions that are not well defined unless they 
are referred to a background spacetime. Standard attempts to introduce a background spacetime tend to use a fixed 
spacetime with properties generally not related to the properties of the system itself, which is commonly considered as 
an unsatisfactory feature (we quite agree with this conclusion). Our procedure, on the other hand, can be considered 
as a universal prescription for defining a unique background spacetime structure driven by, and back reacting on, 
the quantum system. All questions concerning covariance should be referred to this spacetime structure, and this 
spacetime structure is ordinarily not classical. It cannot be considered as a result of a 3 geometry evolution (the 
classical evolution equations do not take a part in our procedure). In particular, we recover the classical evolution only 
through an appropriately-peaked wave function and only through constructive interference over conformal superspace. 

The split of the metric variables into quantized geometrodynamic variables and "classical" or non-quantized variables 
amounts to a representation of a quantum gravitational system as two interacting subsystems. One of the subsystems 
is quantum while the other one is classical. Evolving the many-fingered time variable and coordinatization variables by 
the shift and lapse is just another aspect of such a subdivision of the quantum gravitational system. From this point of 
view the functional evolution problem and the multiple choice problem are simply ill posed. Different time evolutions 
or different quantum theories referred to in the description of these problems are essentially created by introducing 
different physical systems, or by introducing different models of what is believed to be the same system. In other 
words, quantum geometrodynamics yields a finer resolution of different systems than classical geometrodynamics. 
The reason that this split occurs is related to the key features of gravitational systems. In particular, the observer in 
gravitational systems is a part of the system, and the measurement process (coupling of a measuring device with the 
quantum subsystem) cannot be switched on and off at will. In quantum cosmology a particular split of variables into 
quantum and classical essentially reflects om view concerning the nature of time in the Universe. Whether this view 
is correct or not can be decided only on the basis of observations. 

We wish to emphasize here that the square-root Hamiltonian quantization procedure essentially works in a way 
similar to our geometrodynamic quantization procedure. It starts from the split of variables in embedding parameters 
(they are to become the functional parameters of the state functional, essentially classical) and dynamic variables, 
that are to be quantized. A computation of the expectation values of observables involves functional integration 
over the dynamic variables while the embedding parameters are not integrated over. However the approach does not 
contain a meaningful prescription of time parameterization. An attempt to reintroduce such a prescription using shift 
and lapse parameterization of time evolution (cf. Sec. VI) leads to the theory demonstrating the same features as our 
procedure. The spectral analysis problem characteristic to this approach is a consequence of an attempt of one to 
transfer to quantum gravity the fully-constrained property of classical geometrodynamics. Our procedure forfeits this 
feature and instead introduces a system of equations that is not overdctcrmincd, and for which the concept of being 
fully constrained looses its meaning. Nevertheless, this property is recovered in the classical limit only. As a result, 
the geometrodynamic approach does not introduce the square-root Hamiltonian. 

Not every assignment of dynamical variables leads to a consistent quantum dynamical picture. For example, a bad 
choice of the dynamical component of the 3-geometry might lead to a non-elliptic differential operator on the right 
hand side of Schrodinger equation. Furthermore, under such conditions not every shape of the initial wave packet 
will agree with the constraints. In any case, only an appropriate choice of the set of dynamic variables will lead to 
a reasonable quantum mechanical picture. The properties of the ADM procedure together with the results of York 
indicate that there is at least one such reasonable choice. Is there more than one possible quantum dynamical picture 
for the same system? From our point of view, the existing theory is ill equipped to ask this question. The notion of 
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equivalence of quantum gravitational systems is simply not developed. Extreme care should be exercised when one 
attempts to discuss this question naively. 



VI. THE PROBLEMS OF TIME. 



The issues related to time evolution in gravity have been summarized recently by Kuchai'. Q He calls them the 
"problems of time" in quantum gravity. 

The first of the time problems is formulated as a dependence of the final state 'i'spiN ^ given on the 

foliation that connects the spacelike hypersurfaces S/jv and T,pjN. This problem is called the problem of functional 
evolution. 

The second problem is observed by Kuchaf in the context of the square-root Hamiltonian quantization. It stems 
from the fact that the split of the metric variables into the true dynamic variables (3 as opposed to the coordinati- 
zation parameters a and the slicing parameter f2 is not unique. Different splits might lead to different and possibly 
incquivalcnt quantum theories. Kuchai^ calls this the multiple-choice problem. 

The third difficulty is the so called Hilbert space problem and is particular to the Dirac quantization procedure. 
In the square-root Hamiltonian approach, the Schrodinger equation automatically determines an inner product and a 
possibility to form, at given values of a and Q meaningful observables, for which expectation values can be evaluated 
in a way similar to that of standard quantum mechanics (modulo some technical difficulties caused by the functional 
nature of variables). Nothing similar is granted by the Wheeler-DeWitt equation of the Dirac approach. The Wheeler- 
DeWitt equation is a second order functional differential equation for which the problem of supplying the space of its 
solution with the structure of a Hilbert space encounters numerous difficulties. Roughly speaking, the equation does 
not admit a one-Universe interpretation. 

We wish to add to these three problems one more (also described by Kuchaf) , namely the spectral-analysis problem 
that shows up if the square-root Hamiltonian is involved. It is well known that, if the operator expression under the 
square root is not positively definite, the spectral analysis procedure yields, as the result, a Hamiltonian that is not 
self-adjoint. This implies that the Schrodinger equation with such a Hamiltonian does not produce unitary evolution. 

The functional-evolution problem and the multiple-choice problem are, essentially, two aspects of the same problem. 
Both these problems in quantum gravity are caused by (1) the ADM treatment of the whole 3-geometry as dynamic, 
and (2) the fully-constrained feature of the dynamics of the gravitational field. In both procedures of quantization 
(Dirac and ADM) these two requirements are expressed via the demand for the classical constraints to yield operator 
equations that should be satisfied for "any slicing" (or, as they frequently say, for any parameterization). The concept 
of "any parameterization" , however, is essentially classical, which is especially clear when it is thought of as a slicing. 
The concept of slicing works quite well in the classical theory where the 4-geometry of spacetime is unique, but does 
not have any analogue in the quantum theory unless we introduce one and give it an appropriate meaning. We believe 
that both problems emerge as a result of different assumptions concerning the split of a gravitational system into its 
classical and quantum parts which interact with each other. Different splits ordinarily provide physically different 
models of the gravitational system and cannot be compared to each other in the absence of a common spacetime. In 
other words, the concept of equivalent systems is absent in the theory and should be developed prior to posing these 
questions. 

Although there are many different ways to specify coordinatization parameters and slicing parameters, not all of 
them are meaningful in the sense that the resulting time can be related to measurements. A standard procedure to 
introduce a meaningful time in general relativity uses the lapse and shift functions. In this approach the coordinatiza- 
tion and slicing conditions are given arbitrarily only at the initial slice. Subsequently, both the coordinatization and 
slicing conditions are propagated via specifying shift TV* and lapse N functions. In order to facilitate an explanation of 
the consequences of such a description for the square-root Hamiltonian quantization procedure, it is useful to rewrite 
the action principle and Hamilton equations in terms of the /3, a, and Q parameters. 



The expression (2.3) for the action can be written symbolically in the form 

p + n'^a + TT^il - NH{tt'^ ,TT°' ,tt" , (3,a,n) - N,W{tt'^ ,tt" ,TT^ , I3,a,n)] d'^x, (6.1) 



where 



— — -TTT, TT — -, 

S/3 da' dn 



are the momenta conjugate to /3, a, ft respectively. Expression (2.7) takes the form 
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HADAiiT^^, tt", tt", P, a, n, N, N,) = NHiir'^ , tt", tt", /?, a, fl) + N,n'{n'^, tt", tt", /?, a, fi). (6.3) 



The Hamilton equations (2^), ( |2.9| ) become 



(6.4) 



QHadm 



(6.5) 



and 



• OHadm ,„ a\ 



^ = — ' (6-7) 



5/3 ' 



(6.8) 

aa 



■ n dHADM 

TT = 



(6.9) 



The first three equations are merely kinematic relations between the parameters (3, a, Q and their conjugate momenta 
(they can be introduced in any spacetime, not necessarily a solution of Einstein equations). The role played by 



these relations (6.5), is to provide a parameterization (by shift and lapse) of the slicing and coordinatization 

conditions given on the initial hypersurface to the whole spacetime. It is important to realize that this parameterization 
depends on the initial-value problem and is interwoven with the dynamics of the problem. Parameterizing the slicing 
and coordinatization conditions can be achieved by different means. However, one can convince oneself that any 
reasonable set of slicing and coordinatization conditions can be translated into the language of shift and lapse, and 
that such a translation provides a proper interpretation of the conditions in terms of measurements. 



Equations (6.5), (6.6) can be transferred to the square-root Hamiltonian quantization procedure. One should notice, 
however, that the right hand sides of these equations depend on both a, f2 and the true dynamic variables (3 that 
are given, after quantization, by distributions. The situation can be remedied if their expectation values over the 



solution of the Schrodinger equation (augmented by the initial-value data) are substituted in (6.5), (5.6). It is clear 



that no contradictions can be introduced this way (which follows from the simple count of equations and unknowns) . 



The relations (3.5), ( p.q ) can be thought of as definitions of N^, N in quantum gravity. In practice, this procedure 
can become rather complicated. In particular, the parameterization contains functions (the expectation functions of 
dynamic variables as well as of their conjugate momenta) that can be completely determined only after the Schrodinger 
equation has been solved. This means that parameterization of coordinatization functions and the many-fingered time 
scale parameter by the shift and lapse results in an implicit procedure in which even the parameterization itself remains 
undetermined until the final solution of the entire problem has been obtained. Nevertheless, it seems to be unavoidable 
if one is to interpret the solutions of th e quantum gravity equations in terms of observations. Appropriate care should 



be taken, as, in general, equation (6.4) does not have to be satisfied anymore even for the expectations, which means 
that, in general, the expectation values of momenta components are not related to the extrinsic curvature of the 
spacetime determined by the expectation values of the 3- metric, shift, and lapse. The spacetime itself is not a solution 
of Einstein equations. 

A brief summary of the this section is that the problems of time in quantum gravity apparently originate from two 
mutually related sources. 

The first source is a treatment of the classical gravitational field dynamics as dynamics of the 3-metric or 3-geometry 
of a spacelike slice. Both the Dirac and the square-root Hamiltonian quantization procedures are essentially based 
on the original ADM ||^ picture of the gravitational field dynamics. In this picture, the dynamic evolution of the 
gravity field manifests itself as a change from one spacelike 3-geometry to another. In other words, the configuration 
space of geometrodynamics is believed to be Wheeler's superspace. Such an approach does not utilize York's 
analysis of the gravitational field's degrees of freedom and the proper initial-value problem formulation. UJlJ] This 
is not surprising as the foundations of quantum gravity PJlO| were originally formulated before York completed his 
investigation. 
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The second source of difficulties, internally related to the first one, is the demand that the quantum gravitational 
dynamics should be fully constrained, just as the dynamics of classical general relativity. This line of reasoning leads 
to the conclusion that, to quantize the gravity field, one needs only to quantize the constraints. 

The picture of the dynamic evolution of the gravity field as a change from one 3-geometry to another leads 
in classical geometrodynamics to a peculiar situation wherein the Hamiltonian of the Hamilton-Jacobi equation 
(super-Hamiltonian Ti) do not ordinarily coincide with the Hamiltonian participating in Hamilton evolution equations 
CHadm — NH + NiW). It can be argued that other pieces of Hadm are contained in the equations of the super- 
momentum constraint. Nevertheless, the shift and the lapse are thrown out of the picture, and, together with them, a 
meaningful description of a time evolution of gravitational systems. This does not cause any difficulties in the classical 
theory as one can always reintroduce such a description either by retrieving all the Hamilton evolution equations (that 
are compatible with the constraints) or just four of them (the ones that provide kinematic relations between the time 
derivatives of the coordinatization and the many-fingered time functions and the conjugate momenta). One should 
keep in mind, however, that these relations, ordinarily, involve all the conjugate momenta, including those of truly 
dynamic variables. The classical theory on this level treats all the variables in essentially the same way (as functions 
on a spacelike slice). The situation is entirely different for a quantum theory. Some of the variables become the 
arguments of the state functional, others become the functional parameters of the state functional. The arguments of 
the state functional form the superspace. A computation of the expectation value of an observable implies integration 
over the superspace, while the functional parameters are given according to some principles that are not a part of 
the dynamic picture. Ordinarily, the superspace is formed by the dynamic variables. In Dirac's quantization all 
the variables are treated on equal footing and superspace is formed by all the components of the 3-metric. There 
are no functional parameters to be used to fix coordinatization conditions or the many-fingered time. Sometimes, 
as an afterthought, the superspace is reduced to the space of 3-geometries, in which case there are three functional 
parameters to fix coordinatization, but then, in order to satisfy the super-momentum constraints, the observables are 
demanded to be independent of these three functional parameters. The many-fingered time parameter is not in this 
picture. If one decides to treat the many-fingered time parameter as a part of the answer rather than of a question, 
one can try to introduce it as an expectation function of one of the 3-geometry parameters. Questions and difficulties 
related to such a procedure have been outlined in the previous section. 

We wish to emphasize here that this entire situation is created by the treatment of the 3-geometry as the dynamic 
object and an attempt to quantize the 3-geometry. Mathematically, it is expressed via imposing the commutation 
relations on all the components of the 3~metric. One should keep in mind, however, that in general relativity the 
system described by the 3-metric or even the 3-geometry includes in itself the observer and his clock. In standard 
quantum mechanics, or even in the quantum field theory, an observer is external with respect to the quantum system. 
The observer is classical and has an external classical clock. A measuring device can be coupled to a quantum 
system and uncoupled from it at will. There cannot be an external observer in the description of the gravitational 
field because the gravitational system is the Universe itself. An observer cannot switch on and off the coupling of a 
classical measuring device to the system. Dirac's approach essentially quantizes the observer and his clock on equal 
footing with the rest of the system (J. A. Wheeler would say that many-fingered time is quantized). Whether such a 
quantized observer and his clock can function in a fashion providing an opportunity to describe the system consistently 
is not clear to us. A discussion of such a possibility, however exciting, clearly would lead us far beyond the scope of 
this paper. 

Another resolution to the time evolution problem in Dirac's quantization approach would be to first solve the 
Wheeler-DeWitt equation and subsequently split all six variables into (1) the two true dynamic variables /3, and (2) 
the remaining four embedding variables a and f2. The embedding variables could then be reserved for fixing the 
coordinates as well as the many-fingered time parameter, while the true dynamic variables should be integrated over 
only when the expectations of the observables are computed. Unfortunately, such a procedure would contradict the 
entire ideology behind the Dirac quantization procedure. 

The ADM square-root Hamiltonian procedure of quantization avoids the impossibility of introducing a functional 
time via solving the super-momentum and super-Hamiltonian constraints on the classical level and quantizing the 
resulting equations. The superspace of the theory is thus reduced to the space of true dynamic variables with 
four functional parameters left as coordinatization parameters and a many-fingered time parameter. This procedure, 
however, leads to a Schrodinger equation with a square-root Hamiltonian, which poses well known technical difficulties 
briefly discussed above. That no prescription for a meaningful choice of the coordinatization and many-fingered time 
variables, similar to the shift and lapse parameterization of classical general relativity, is supplied in this approach 
presents another difficulty of a conceptual nature. An attempt to augment the procedure by such a prescription is 
described above. The description of the evolution becomes implicit with a parameterization being finalized only after 
the entire problem of evolution has been solved. We wish to stress at this point that the ADM square-root approach to 
gravity quantization starts, as Dirac's approach, from a treatment of classical gravity field dynamics as the dynamics 
of a 3-geometry. Exactly as in Dirac's approach, the fully-constrained property of classical general relativity is 
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introduced as a key feature that is to be transferred to the quantum theory Uterally. This is accomphshed by turning 
the constraint equations into operator equations which restrict the admissible states of the gravitational systems. The 
only difference is that the constraints are solved before quantization in order to introduce a Hamiltonian that, in a 
sense, can be interpreted as a time evolution generator. The (functional) dimension of the superspace is reduced to 
the correct dimension as determined by a proper analysis of gravitational degrees of freedom. [^Q] However, this 
achievement comes almost as an afterthought, a fix of a deficiency that has been created by the interpretation of the 
3-geometry as the principal dynamic object, which, in turn, has forced on the theory an identification of the constraints 
with the Hamilton- Jacobi equation. The square-root Hamiltonian participating in this equation again differs from 



the Hamiltonian Hadm participating in Hamilton equations (2.8), ( |2.9| ) and is not related to the Lagrangian by a 
relation similar to (^.15 ). Such a switch of Hamiltonians is possible only because the classical dynamics of general 
relativity is fully constrained. In transition to the quantum theory, the achievements of this entire procedure come at 
a price, creating difficulties on both the technical and the conceptual levels. 



VII. DISCUSSION. 



In this paper we have reviewed the issue of time evolution in quantum gravity. Both the Dirac and the ADM square- 
root Hamiltonian quantization procedures create difficulties in introducing a meaningful concept of time evolution. 
The difficulties in both approaches, we conclude, stem from their common tendency to transfer to quantum theory 
an interpretation of classical geometrodynamics as an evolution of a spacelike hypersurface 3-geometry together 
with the fully-constrained property of classical geometrodynamics. These two features essentially lead in classical 
geometrodynamics to an identification of Hamilton- Jacobi equation with the constraints and thereby removing a 
meaningful concept of time. This procedure does not present a serious problem in the classical theory of gravity as 
time evolution can be reinserted back. However, it leads to two different Hamiltonians in the theory. One of them 
participates in Hamilton evolution equation, the other one is the Hamiltonian of Hamilton- Jacobi equation. 

These two features of classical general relativity, when transferred to quantum gravity, motivate quantization of 
constraints (in original form in Dirac's procedure, or resolved with respect to the momenta conjugate to the embedding 
variables in the ADM square-root Hamiltonian procedure) via imposing their operator versions on the state functional. 
The resulting quantum theories naturally counter any meaningful concept of time evolution. They lead to what Kuchai' 
has identified as the "problems of time," including (1) the problem of functional evolution, (2) the problem of multiple 
choice, (3) the Hilbert-space problem of Dirac's approach and (4) the spectral-analysis problem of the square-root 
Hamiltonian approach. 

We have discussed the problems of time mainly for the ADM square-root Hamiltonian quantization, as Dirac's 
quantization, when interpreted literally does not admit a functional many-fingered time. Presumably, it can be 
somehow reinterpreted (cf. Sec. II), but one such avenue leads back to the ADM square-root Hamiltonian approach. 
This approach has all the time problems in it except the Hilbert-space problem. The way the ADM approach is 
ordinarily formulated leaves little room for an analysis of the functional-evolution problem and the multiple-choice 
problem. At first sight, these problems appear to be genuine. The spectral-analysis problem seems to be unavoidable 
for the square-root Hamiltonian quantization. In addition, we conclude that the approach suffers from the lack of a 
meaningful prescription for time parameterization (although it admits functional time). Our attempt to reintroduce 
into the ADM quantization a time parameterization of evolution by shift and lapse has lead to an implicit set of 
equations with a parameterization depending on the initial values. This last feature does not destroy the procedure 
but it gives one an idea of what one is to expect if he introduces time in a way accessible to measurements. 

The new idea presented in this paper is to reconsider the classical picture of geometrodynamics via representing 
it as a dynamic theory of two interacting subsystems, one being described by the dynamic components of the 3- 
metric, and the other being described by the many-fingered time variable and coordinatization variables. If one is 
to keep the concept of classical time in the theory, one should quantize only the first subsystem, while treating the 
second one as classical yet generated by, and back reacting on, the first (quantum) subsystem. In such an approach 
classical geometrodynamics can be described by two sets of equations. Equations describing the evolution of the 
first subsystem treat the second subsystem as an external field. Having in mind that, eventually, we are to quantize 
the first subsystem, we describe its evolution by an appropriately formulated Hamilton- Jacobi equation. For the 
second system we use as field equations super-momentum and super-Hamiltonian constraints rewritten in appropriate 
variables. The full system of equations is equivalent to the ten Einstein equations. However all the equations of this 
system are independent. 

Quantization of such a system is achieved via turning the Hamilton-Jacobi equation into a Schrodinger equation. 
The constraint equations describe the second subsystem as classical but include the dynamic variables as a source. 
We replace them by the expectation values of these variables over the solution of our Schrodinger equation (with 
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appropriate initial data). It is clear that in such a procedure of quantization there is no Hilbert space problem. The 
Hamiltonian of the Schrodinger equation does not contain square roots which eliminates the spectral-analysis problem 
as well. 

The problem of functional evolution and the multiple-choice problem obtain an interesting interpretation in our 
proposed approach. In our gcomctrodynamic procedure of quantization they turn into a statement that the solution 
of a given problem depends on the split of the original system into its quantum and classical subsystems. In par- 
ticular, different splits generate essentially different systems. In other words, quantum geometrodynamics resolves 
gravitational systems finer than classical geometrodynamics. 

We have illustrated our procedure using a Bianchi lA cosmology and an axisymmetric Taub model. When these 
models were quantized, our theory generated quantum geometrodynamic pictures based on a post-ADM treatment 
of the gravitational field dynamics and was free from the conceptual difficulties usually associated with the Dirac 
and ADM procedure of quantization. The variables describing the gravity field in these cases have been split into 
the true dynamic variables and a parameter related to Wheeler's many-fingered time. Only the dynamical part, 
the underlying anisotropy, has been quantized. The Hamiltonian participating in the Schrodinger equation is not a 
square-root Hamiltonian. This absence of a square-root Hamiltonian is generic for our quantum geometrodynamic 
procedure. The effective "background spacetime" determined by the expectation values of dynamic variables together 
with the "observers" (related to Wheeler's many-fingered time) allows us to pose unambiguously the questions of 
covariance. 

The problems outlined in Sec.Vl become all but eliminated by our quantum geometrodynamic approach. The 
nontrivial part of gravity quantization appears to shift from such conceptual problems too the problem of (1) the 
choice of an appropriate model to quantize, and to the related problem of (2) the choice of an appropriate initial 
condition for the wave functional. Both choices are crucial if one is to attempt using quantum geometrodynamics 
to better comprehend the properties of gravitational systems. It is our understanding that the success of gravity 
quantization rests on such meaningful choices. Furthermore, the choice of models should not be determined by the 
structure of quantum geometrodynamics; rather, it should be determined by observational data and our general 
understanding of gravitational phenomena. 

It is clear that the procedure of quantizing the dynamical part of constrained systems described in this paper can be 
performed in the general case of geometrodynamics without any complications in principle, although it may become 
quite involved computationally as compared to our simple model examples. p8[ | The procedure differs only in two 
respects from the simplistic examples presented here. The first difference arises when the 3-metric is (1) parameterized 
by three coordinatization parameters, (2) the many-fingered time parameter, and (3) the two dynamic variables; then 
all four constraints should be solved with respect to the coordinatization parameters and the scale factor. In all four 
constraints the expectation values of the true dynamic variables should be used. The second difference is caused by 
the functional nature of the gravitational field dynamics in the general case. The operation of functional integration 
is involved, which might lead to analytic difficulties. Such difficulties are not specific for our approach as they are 
common for the canonical formulations of all field theories. Quantum geometrodynamics, in particular, does not seem 
to generate any specific new difficulties. In this paper we forwarded the beginnings of a quantization scheme consistent 
with York's analysis of the gravitational degrees of freedom. Our particular imposition of the four constraint equations 
leads to a weaker theory that in turn avoids the problems of time. 
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